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Abstract. A finite (periodic) FPU chain is chosen as a convenient point for 
investigating the energy exchange phenomenon in nonlinear oscillatory systems. As 
we have recently shown, this phenomenon may occur as a consequence of the 
resonant interaction between high-frequency nonlinear normal modes. This interaction 
determines both the complete energy exchange between different parts of the chain and 
the transition to energy localization in an excited group of particles. In the paper, we 
demonstrate that this mechanism can exist in realistic (asymmetric) models of atomic 
or molecular oscillatory chains. Also, we study the resonant interaction of conjugated 
nonlinear normal modes and prove a possibility of linearization of the equations of 
motion. 

The theoretical constructions developed in this paper are based on the concepts 
of "effective particles" and Limiting Phase Trajectories. In particular, an analytical 
description of energy exchange between the "effective particles" in the terms of non- 
smooth functions is presented. The analytical results are confirmed with numerical 
simulations. 

Keywords: Nonlinear oscillatory chain; Nonlinear Normal Modes; energy transfer; 
localization; "effective particles"; Limiting Phase Trajectories 
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1. Introduction 

It is well known, that the problem of vibration energy localization in the infinite 
Fermi-Pasta-Ulam (FPU) chain can be asymptotically reduced to finding a localized 
solution (breather) of the continuum Nonlinear Schrodinger Equation (NLSE), which is 
a completely integrable system [lj. After such reduction to NLSE, the mathematical 
origin of the phenomenon becomes clear. However, it was recently shown that not an 
infinite but a finite (periodic) /3-FPU chain has to be chosen to clarify the physical 
nature of vibration energy localization [2]. This model has been intensively studied 
over the last decade ([3], [3], [5], [6], [7]) in the framework of Nonlinear Normal Modes 
(NNMs) concept [8], [9]. On the contrary to NLSE, finite discrete systems are not 
completely integrable (except the three particle /3-FPU periodic model, [in])- Therefore, 
the knowledge of NNMs leads to a unique possibility for constructing a wide class 
of solutions using a perturbation techniques pp. We have recently shown that an 
increase of the particles number in the symmetric ( /3-FPU) chain leads inevitably 
to the resonant interaction between the high frequency NNMs because of the linear 
spectrum densification. We suggested the concepts of "effective particles" and Limiting 
Phase Trajectories (LPTs) to understand and describe both complete energy exchange 
between different parts of the chain and transition to energy localization [2]. From 
the modal viewpoint, such exchange results from the resonant interaction of the zone- 
boundary mode (7r-mode) with its nearest neighbours. In this case, the formation of a 
mobile localized mode is accompanied by a prohibition of the transfer of energy from 
one part of the system to another. The complete energy exchange between "effective 
particles", which is described using LPT concept [13], is an alternative to nonlinear 
normal oscillations characterized by the conservation of energy. 

The symmetric model of the FPU chain is valid for various mechanical problems. 
However, for applications to the problems of Molecular and Polymer Physics it is 
necessary to consider the asymmetric (a — /3) FPU problem because realistic interatomic 
potentials of interaction are strongly asymmetric. The presence of the asymmetric term 
in the interaction potential leads to significant complication of the problem. First of 
all, the NNMs of the discrete FPU chain with periodic boundary conditions cannot be 
divided into closed 4- manifolds, as it is possible in the case of the /3-FPU chain [2] . The 
modal equations for the degenerate pairs of modes corresponding to the same eigenvalue 
are valid only for special initial conditions. 

One of the main points under consideration is whether the presence of the quadratic 
terms in the equations of motion leads to an effective renormalization of the constant of 
the fourth-power interaction. Obtaining the positive answer, we present an analytical 
description of energy exchange between "effective particles" in the (a — /3)-FPU chain 
using the LPT concept. 

Note that the same effective renormalization of the quartic coupling constant was 
introduced for modulational instability analyses of modes close to the zone-boundary 
mode in the rotating- wave approximation ([H], [12]). The current paper considers a 
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more wide class of problems, a the key result of this paper is in an analytical description 
of intensive energy exchange between "the effective particles" and energy localization 
on the effective particle in terms of the non-smooth basic functions. 

We also study the asymmetry effect on the threshold to energy localization. 
In conclusion, the comparison of analytical solutions with the results of numerical 
simulation is performed. 

2. The model and the equations of motion. 

We consider the (a — /3)-FPU system defined by the Hamilton function 

N I I a p 

H ° = Y1 2^ + 2 ~ 9j ^ + 3 ( Qj+1 ~ + 4 ~ ^ 

3=1 

and the periodic boundary conditions (qN+i = Qi,Pn+i = Pi), where qj and pj are the 
coordinates and the conjugate momenta, respectively, N is the number of particles. 

The transformation to normal coordinates is given by the linear canonical 
transformation 

N-l 
fc=0 

with the coefficients 

^ = _ (an — + cos — j , k = 0,l,...,N-l. 

The transformation (2) allows us to present the quadratic part of Hamilton function 
as the total energy of the independent oscillators: 

N-l 

H2 = J2^t + cote k )- (3) 

Here ^ and r\ k are the amplitudes and the momenta of NNMs; the coordinate £ 
associated with motion of the center of mass is removed from (3). The eigenvalues u k 
have the form 

Wfe = wat/2 sin— , wat/2 = 2, k = 0, . . . , TV - 1. (4) 

If the number of particles increases, the frequency gap between the highest- 
frequency mode and nearby modes quickly decreases. In what follows, we consider 
an even number of particles N. In this case, the frequencies u k are bounded by the 
highest-frequency oujy/2 = 2. The zero eigenvalue uq = corresponds to motion of the 
chain as a rigid body. All eigenvalues (4) are in the interval < u>k < 2 and twice 
degenerate. 
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The remaining part of the Hamilton function (1) for the periodic (a — /3)-FPU chain 
has the following form [3]: 

N-l 

a 
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i,j,k=l 

N-l 



^k^l^m^nik^lim^nCklmn^ (5) 



where 



and 
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otherwise 



So the equations of motion can be written as: 



,2C N - 1 



d £k 2 C a 

—rrr + W fc ^fc H ~F^ u k / WiU m Di m kqiq m + 

at 1 2ViV 

l,m=l 

Jr ^2N° Jk ^l^m^nCklmnilimin = 0. (6) 

l,m,n=l 

To analyse the dynamics of the chain in the context of Eqs (6), we introduce the 
complex amplitude corresponding to the combination of displacement and velocity in 
the new basis 



1 



( a* 



^k = -j= I ^ + iX k£kJ , u k =: X 2 k , 
In terms of the complex amplitudes, Eqs (6) can be rewritten as 

, T N-l 

iu k ^! k = \ X k Xi\ m D lmk {% - V*) (* m - V* m ) + 

A * /O AT < * 



dt 4V2iV , , 

l,m=l 

+ |v ^ ^kXlX m X n C k imn (*i ~ (*m ~ *m) (*n ~ *«) = 0. (7) 

l,m,n=l 

Taking into account the dependence of the chain properties on the number of the 
particles, one can choose the value l/\f~N as a small parameter e. 
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3. Main asymptotic approximations. 

We employ the multiple scale procedure [8] in order to consider the processes in the time 
intervals greatly exceeding 2-n/uk. Following this procedure, we introduce time scales: 

T = t, Tl = £T , T 2 = E 2 T , ... (8) 

where the "fast" time r corresponds to the initial time scale of the system, while the slow 
times ri, r 2 , etc. correspond to the slowly varying envelopes. The envelope functions 
(fk are defined by the relation 

We construct the asymptotic representation of tp^ in the form 

= Xk,i + £Xk,2 + £ 2 Xk,3 + ■■■ (9) 

The multiple scale expansion based on the relations (8), (9) will be used to derive 
the equations of the leading-order approximation Xk,i- 

After some calculations one can get the equations for the pair of modes with the 
same eigenvalue (we denote Xk,i ='■ Xk and suppose that the other modes are not 
excited.): 

U + c ltk \xk\ 2 Xk + c 2 , k \XN-k\ 2 Xk + c^ k x 2 N _ k xl = 



dr 2 

' <h "^~ k + c 1>k \xN-k\ 2 XN-k + C 2 ,k\Xk\ 2 XN-k + c^ k x 2 k X*N-k = °) ( 10 ) 



dr 2 



where 



9(3 3a 2 \ 2 2 \ //3/3 a 2 



c hk = [ I ^ — ) + " 2 I , °2,k = ( [ — + — )u} 2 k -2a 2 



4 2 



3/3 5a 



|2 , |,, |2 



C3,k = ^ ^ — — ) + 3a 2 

This system possesses integrals 

X = \xN-k\' z + \Xk\ 
and 

G- = i(XkX*N-k ~ XN-kXl), 

so it can be rewritten as: 

dXk 



dr 2 
dxN-k 



+ c 3tk G-XN-k - ic lik Xxk = 
C3,kG-Xk - ici^XxN-k = 



dr 2 

The effective linearization of the equation of motion (10) clarifies why the normal 
modes may preserve their individuality in the nonlinear case. The intensive energy 
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exchange corresponds to small value of G_. This process is described by a phase 
trajectory which is far from the equilibrium point, and results in a slow rearrangement 
of energy by periodic waves. The rate of energy exchange increases as the value of G_ 
grows, but amount of the energy participating in the exchange decreases (Fig.l). 



X = 1; G_ = 0.61 




X = 1; G_ = 0.79 



X = 1; G_ = 0.97 



£0.06 
0.0+ 




2 0.06 
0.0+ 
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Fig.l. The increase in energy exchange rate and the decrease of the energy exchange 

intencity while growing GL. 

The dynamics of the normal modes in the essentially discrete system (small value 
of N) should be considered in the assumption that the gap between 7r-mode and nearby 
conjugate ones is large enough. But, if the number of particles grows, then the new 
opportunity arises. Namely, the resonant interaction between the modes having close 
frequencies (in the linear approximation limit) becomes possible. 

While the energy of excitation and/or number of particles grows, the dynamics of 
the chain becomes nonquasiharmonic. While the parameter e decreases, the eigenvalues 
of the 7r-mode and the nearby ones become more and more close until the resonant 
conditions occur. 



W N/2-1 — ^N/2 



1 - 



71' 



2N 2 



In such a case the equations (10) are not valid because this resonance is not taken 
into account. Now we choose the top frequency un/2 as the basic one. So, because of 
the closeness of frequencies, additional linear terms appear in the equations for Xn/2-i 
and Xn/2+i- Moreover, we have to take into account the resonant interactions in the 
nonlinear part of equations. The problem turns out to be more complicated than that 
for the /3-FPU chain. After lengthy transformations (details of calculation one can see 
in Appendix) we have got the the equations 



■ gXi 
dr 2 

; 9X2 

' dr 2 



k r>2 



7T 



2 

■ dxo 



tf X 2 



+ 2CVL 



Cn 2 (\ X i\ 2 Xi + xlxl + 2xlxl + 2\xo\ 2 Xi) = 

CQ 2 {\x2W2 + xlx*2 + 2xlx*2 + 2\xo\ 2 X2) = 
2 ((\Xi\ 2 + \X2\ 2 )Xo + (xl + xl)x*o))=0, 



11) 



where 



C 



3/3i 
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3/3 
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Xo 



Xn/2, Xi '■— Xn/2-i, X2 '■— Xn/2+i, ^ ^>N/2, 
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which coinside with those for the /3-FPU chain if the constant of the nonlinear interaction 
is renormalized (/3 fa). This renormalization is in a good agreement with the 
numerical results, e.g. the asymptotic system in slow time obtained from FPU chain with 
such a and fa that fa — behaves almost like linear one (Fig. 2). So all effects described 
in [2] also occur in a/3-FPU chain. Therefore, in the framework of this asymptotic 
approximation there is no difference between a fa and /3-FPU chains (the dynamics of 
a-FPU system corresponding to < 0). But, contrary to [2], we do not restrict ourselves 
by numerical results and perform the analytical solution of (11) for LPT. Besides, we 
discuss in more details the notion of the "effective particle". 




Fig. 2 Distribution of the energy among the particles of the chain; initial condition 
corresponds to LPT (excited one half of the chain - "effective oscillator"): a) linear 
system (a = (3 = 0) and b) "effectively linear" system (a = 0.5, = 0.333). 

It is easy to check that equations possess an integral of motion - total "occupation 
number" : 

Xtot ■= \xi\ 2 + |X2 1 2 + |xo| 2 - 



4. From "waves" to "particles". 



The description of system's dynamics in the terms of NNMs has to be replaced by that 
in terms of the "effective particles". This allows us to describe the intensive energy 
exchange between different parts of the chain, which can be identified as "effective 
particles" with "coordinates" ipi and ip2' 



4>i 



1 

7! 



Xi 



Vl - C 2 + CX2 



<* 2 = 71 



Xo + (xiVl - C 2 + C\2 



V 3 = [CXl ~ X2 1 

where c is a constant defined by initial conditions (0 < c < 1). Such transformation 
preserves the total occupation number in the form X = \tpi\ 2 -\-\ip2\ 2 +\v\ 2 - For the value 
c = which corresponds to equal initial conditions of modes, we have the equations 
in the coordinates of "effective particles": 



Mi , 7T 2 2 , 

l W 2 + T fi ^ 



Vi) + CQ 2 \i> 



9 / 
4^ 



1^2 



3 , 1 , 

2^1 + 4^2 
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Since the notion of effective particle is of great importance, we discuss it in more 
details. When dealing with the mechanisms of intensive energy exchange, one needs 
to reveal the elementary agents which exchange energy and become the domains of its 
localization (after exceeding some critical level of excitation). In gaseous media, they are 
the interacting particles (atoms or molecules) themselves, which participate in almost 
free motion. 

On the contrary, the particles in oscillatory chains as well as in all crystalline 
solids undergo a strong mutual interaction. Therefore, the elementary agent here is all 
oscillatory chain performing oscillations corresponding to one of non-interacting normal 
modes. However, an increase of the particles number leads to the resonance relations 
between several Normal Modes and then between multiple modes with close frequencies. 
If the initial conditions are strongly asymmetric (they are far from those corresponding 
to every of resonating normal modes), it leads to the manifestation of two significant 
effects. 

First of them is the inter-modal coherence and second one is the inter-modal 
interaction (in nonlinear case). As a result, the resonating normal modes are no longer 
the appropriate elementary agents. In the system of two weakly coupled oscillators 
their role is played by the particles themselves (as well as in the case of gaseous 
media) changing slowly by the energy in the beating process (their displacements can 
be presented as a sum and a difference of the modal variables). 

As the number of particles increase, the appearance of the resonating modes does 
not mean that we have to come back to the real particles. Now the "effective particles" 
play a role of elementary agents. 

Their displacements can be constructed as combinations of resonating modal 
variables only (similarly to the displacements of two weakly coupled oscillators) which 
manifest the inter-modal coherence in the beating process described by limiting phase 
trajectories. Thereby, the beating notion is extended to multidimensional systems. 
Besides, the introducing the LPT allows us to describe adequately the transition from 
intensive energy exchange to energy localization due to an increase of intensity of 
excitation. By this manner the connection with continuum systems, having the breathers 
as localized solutions is clarified. At last, a most simple analytic presentation of intensive 




+ 




(12) 
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energy exchange can be attained precisely in terms of the "effective particles" and 
Limiting Phase Trajectories. 

Let us consider a particular solution corresponding to ip — 0, xi — Xi- I n this case, 
the Hamilton function is 

H = Q 2 (-^1 - H 2 + C (9(|^| 2 + |^ 2 | 2 ) 2 + Q\H 2 \H 2 - 

-2(i^i 2 + \H 2 )(^r 2 + mi) + tfm 2 + itm 2 ) ) . as) 

Since the occupation number X = \ipi\ 2 +\ip 2 \ 2 i s the integral of motion, the variables 
■01 and ip2 can be expressed through the angle variables 9, 61,62: 

ip 1 = X cos 9e iSl , ip 2 = X sin 9e i&2 (14) 

Substitution of (14) with ip = into (12) yields the following equations 
d9 

h Ki sin A + K 2 sin 29 sin 2 A = 

dr 2 

^ sin 29 + 2K X cos A cos 29 - AK 2 sin 29 cos 29 (8 - cos 2 A) = (15) 
dr 2 



where 



A = 6 1 -6 2 , K x = n 2 ^— ; X 2 = fi 2 ; ft = p - -a 2 = -C. 



4 32 J 1 64 ' " ' 3 3 

Note that Eqs (15) conserve the first integral 

if (0, A) = X ( 27/3lX 6 ~ 16n2 ^ 2 + K, sin 20 cos A- 

-K 2 (8 - cos 2 A) sin 2 20) (16) 

obtained from (13) by changing ipi,ip2 to expressions (14). 

Equations (16) for all X, ft have the following stationary points (Fig. 3): 

(a) A = 0, 9 = £ (6) A = tt, = J (17) 

The points (a) and (6) correspond to ipi = ^ 2 ', by definition, this leads either to the 
pure 7r-mode (a) or to the pure (xn/2-i + Xiv/2+i))-mode (6). 
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Fig. 3 Phase portrait of Eqns (15) with X < X c (LPT - dashed line) . 

When the amplitude of excitation grows, the in-phase stationary point (17a) 
becomes unstable (Fig. 4). This instability leads to the appearance of two new stationary 

2 

points. Namely, if the occupation number X exceeds the value X c = then there exist 
three stationary solutions of Eqns (15) with the phase shift A = 0: 




Fig.4 Distribution of the energy among the particles of the chain (numerical solution); 
initial condition corresponds to 7r-mode: a) before the treshold of instability X c and b) 

after it. 

The threshold of instability of the 7r-mode is equal to X c = 7r 2 /3/?i, which is in a 
good accordance with the estimation obtained in the framework of the "narrow packet" 
approximation [15J. Any trajectories in the neighborhood of stationary points (186, 18c) 
correspond to weakly localized solutions, for which the energy of one half of the chain 
- "effective particle" - only slightly exceeds the energy of the second one. At the 
same time, the trajectories starting beyond this vicinity and passing through the points 
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ipi (9 = 0) or ip 2 (9 — tt/2) correspond to the combinations of modes with approximately 
equal energies that means that an initial excitation of one "effective particle" entails 
the complete energy exchange between both ones, i.e., the transition from the state 
ipi to the state ip2 and inversely. This implies that a possibility of complete energy 
exchange between different parts of the chain exists for the excitation level exceeding 
the instability threshold X c (Fig. 5). 




Fig. 5 a) Phase portrait of Eqs (15) (LPT - dashed line) and b) distribution of energy 
among the particles in the chain along LPT (numerical solution), X c < X < Xi oc , . 

A growth of the amplitude X entails an enlargement of the domain encircled by the 
separatrix passing through an unstable stationary point (out-of-phase NNM); at last, 
the separatrix coincides with the LPT. At this instant, the topology of the phase plane 
changes drastically (Fig. 5). Namely, any energy exchange between the mixed states 
ipi and ip2 disappears; this implies that a trajectory starting at a point corresponding 
to 9 < 7r/4 (or 9 > vr/4) and for any A cannot reach a point corresponding to 9 > it/ 4 
(or 9 < 7r/4) [excepting the trajectories surrounding in-phase stationary points 9 = ±7r 
and bounded by the separatrix going through the unstable point 9 < 7r/4, A = 0]. 
Therefore, the energy initially concentrated near the states ip\ or ip2 remains confined 
in the excited effective oscillator (Fig. 6). 




Fig.6 a) Phase portrait of Eqns (16) (LPT - dashed line) and b) distribution of energy 
among the particles in the chain along LPT (numerical solution), X > Xi oc . 

An energy threshold associated with the above mentioned localization can be found 
from the condition of equality of the energy corresponding to the LPT and the energy 
at the unstable 7r-mode. It is seen in Figs. 3, 5, 6 that the LPT goes through the points 
9 = 7r/2 and 9 = 0. This means that 

H(9, A)\ LPT = ^ - 16tt 2 ) 
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cos A = ^ sin 29 (8 - cos 2 A) , 



4 32 y 64 

Since the energy of the 7r-mode is equal to zero, it is easy to calculate the respective 
occupation number X; oc = 167r 2 /27/3i and the energy of the chain E[ oc = 167r 2 /27/3iiV. 
It now follows that above the excitation level E ioc we can observe the localized vibration 
excitation (a breather). 

If the initial energy is concentrated at the state ipi or ifj 2 , the representing point 
in the phase plane moves along a trajectory encircling the respective stationary point. 
Then the temporal evolution of the breather corresponds to the regular variation of its 
profile (the "breathing" mode of the localized excitation). The period of breathing can 
be calculated as 

where the last integral is taking along the LPT. 

Contrary to "breathing" breathers corresponding to the motion along the LPT, the 
stationary points (19) determine new normal modes with invariable non- homogeneous 
energy distribution, or the breather-like excitations. In spite of these modes existence 
at any X > X c , the true threshold of localization is equal to X/ oc because the possibility 
of the complete energy exchange is preserved for motion along the LPT till X = Xi oc . 

Thus we get the solution (18) corresponding to immobile breather. Finally, we 
want to clarify the nature of travelling breathers. We recall that the above results 
have been obtained under assumption ip — that corresponds to equal amplitudes of 
the conjugate modes. Now we assume that the amplitudes of the conjugate modes are 
slightly different, so we consider small enough ip ^ 0. One can show that the equations 
for functions ip\ and ip2 include only quadratic terms depending on (p. Therefore, in the 
framework of the linear approximation, if the value <p ^ is small enough, there is no 
qualitatively change in the phase-plane portrait in Figs (3, 5a, 6a). As for the behaviour 
of ip in the vicinity of any of the stationary points (18), it is described as follows: 



- % + ^M( 2 ^sin2# + (3 + e 2 ^sin26)p*) = 0. 

GST2 2 8 



The respective eigenvalue 



v/891/3 2 X 2 + 607r 2 /3!X - 96vr 4 

A — it 

28 



is imaginary if X < Xi oc , i.e. if the amplitude of excitation is less than the threshold 
corresponding to coincidence of LPT and separatrix. 
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Fig. 7 Directional motion of the breather when a real part in the eigenvalue A at 

X > Xi oc . 

It is easy to demonstrate that small non-zero ip in the equations of motion leads 
to small oscillations of the breather center when the eigenvalue A is imaginary. From 
the other hand, appearance of a real part in the eigenvalue A at X > X\ oc leads to the 
directional motion of the breather (Fig. 7). 

5. Analytical solution for the LPTs. 

Using the equation for LPT 

K 

cos A = — sin 29 (8 - cos 2 A) , 



one can get the solution for LPT in terms of elliptical integrals: 

— - 2yfK$ (8 - cos 2 A) 2 - K\ cos 2 A = 0, 

A 

dA, 

r 2 



2 

cos 2 A 



r 2 



2\jKl (8 - cos 2 A,) 2 - K\ cos 2 A, 
dz 

VO? - z 2 )(Kl(8 - zf - Kffi' 



We get the formal solution of (15) to LPT, but, since its use is prohibitively 
difficult, we analyse the system in angle variables (15) with LPT equation with 
X < Xi oc , X > Xi oc and derive an asymptotic formula to a periodic solution of these 
equations. 
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It is seen in Fig. 8 that 9{j 2 ) (X < X [oc ) and A(r 2 ) (X > X ioc ) forms are close to the 
straight line, with an almost instant reverse at n/2 in the case of 9 and instant step at 
— 7r/2 in the case of A. In order to describe the solution we introduce the new variable 
t(<P), <p = Q*t 2 , - constant, which will be yielded from the first order approximation. 
The function r(0) and its derivative e(0) = dr/d<p have the form 

r(0) = — | arcsin(sin — )|, 

7T 2 

e(<f>) = 1, < <p < 1; e(» = -1, 1 < <p < 2. (19) 

Following [9], we construct an approximate solution as a function of r; the inverse 
transformation r 2 = t 2 (t) automatically yields a periodic solution in r 2 . Taking into 
account the discontinuity of A(r 2 ), we construct the solution of (15) in the form 

0(T 2 )=e(r), A(r 2 ) = e (0)F(r), A = n, e iL (20) 

ar 2 ar 

To derive the equations for 0(t), V(t), we insert (19) into (15). This yields the set 
of equations 

<96 

ft,— + K lS mY + ^2 sin 26 sin 2F = (21) 
or 

QY 

fi* — sin 26 + 2Ki cos Y cos 29 - 4K 2 sin 26 cos 26(8 - cos 2 Y) = 
or 

e = e(r(sj»T 2 )), e(o) = o; y = r(T(n,T,)), y(o) = ~, 

A - 1= tf(^-^;/ f2 = ^tf. 
Using the equation for LPT 



4 32 J 1 64 



cos F = ^ sin 26 (8 - cos 2 Y) , 

one can get the leading-order approximation of the solution of system (21) in the case 
of beating between "effective particles" far from localization: 



6 (r) = A r, A + Ax = -| + arctan f ^-Htan ( sin2.4,,,r j | . [22) 



o"*o 

A = ^,A n, = K x 
and in case of strongly localized solution: 

6 = A* sin 7tt, y = — | r |,A*= arcsin — — A^fi^n = Xi . 

2 V 2 J 7K 2 

Fig. 8 demonstrates that the numerically constructed solution of (18) is in a good 
agreement with the yielded approximations. 
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Fig. 8 Comparing of the numerical solution of (15) with the analytical first order 
approximations in the case of beating between the "effective particles" (a, b) and in 
the case of strongly localized solution (c, d). 

6. Conclusion 

This study demonstrates that the energy localization phenomenon in the nonlinear 
oscillatory chains can be described in the framework of the unified concept for the 
systems with both symmetric and asymmetric interaction potentials. This concept is 
based on the notions of "effective particle" and Limiting Phase Trajectory. It has 
been shown that asymmetric interaction, typical for molecular systems does not change 
qualitatively the main features of the energy localization phenomenon observed in the 
symmetric nonlinear oscillatory chains. In fact, the revealed renormalization of the 
four-order interaction constant substantially simplifies the further analysis. Besides, 
this allows us to conclude that asymmetry of the potential leads to reducing the 
effect of nonlinearity on the dynamical properties of the system. However, when 
this renormalized constant becomes small, the higher-orders terms in the expansion 
of the exact potential energy have to be taken into account. The transition to energy 
localization turns out to be a consequence of the prohibition of the intensive energy 
exchange between "effective particles". The latter is an adequate analogue to the 
beating phenomenon, well-studied in the systems of two oscillators. The analytical 
results for Limiting Phase Trajectories describing the energy exchange and localization 
are confirmed with numerical calculations. 

Appendix A. Multiple scale procedure. 

We employ the multiple scale procedure [14] in order to consider the processes in the 
time intervals greatly exceeding 2tc /u^- Following this procedure, we introduce the time 
scales: 

r = t, n = et , t 2 = e 2 r , . . . (Al) 
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where the "fast" time r corresponds to the initial time scale of the system, while the slow 
times Ti, r 2 , etc. correspond to the slowly varying envelopes. The envelope functions 
ipk are defined by the relation 

di Pk clulk t \,\.\ n. . _ 



dt 4 



Lm=l 



Bi N ~ 1 

i,m,n=l 

* (</? m e^ - *) (^ n e^* - ^e"^*) = 0. 

We construct the asymptotic representation of ip k in the form 

Vk = Xk,i + £Xk,2 + £ 2 Xk,3 + ■■■ (A2) 

The multiple scale expansion based on the relations (Al), (A2) will be used to 
derive the equations of the leading-order approximation Xk,i- 
Terms with e° : 

d Xk,i = Q 
dr 

Terms with e 1 : 

N-l 



^ 1 ' Lm=l 



* (x™,ie™ - x^e"™) = 

There is one secular terms in this equation: = =>• Xfc,i — Xfc,i( r 2)- Then we can 
get the expression for Xk,2 and use it at the further steps of approximation: 

e i(u> p -ui q -u; k )T g-i(wp+w q +w fc )ro \ 
- 2 Xf>,lXq,l 7 r - Xp,lXg,lT7 ; ; r • (A3) 

Terms with e 2 : 

e iu k r ( d Xk,s + d Xk ,2 + <9Xfc,l 



<9r 9ri <9t 2 



«V2 ^ 



x ^k\r\sD rsk {{ X r,2Xs,l + Xr,lXs,2)e^ + ^ ~ 

r,s=l 

-l{x* r ,iXs,2 + xhXs,i)e^'^ + (x^xh + X*r,2X:,i)e- i{ ^ )T0 ) + 
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■ a N-l 
l,m,n=l 

* {Xm,ie^ - x*m^ mT0 ) (Xn,ie™ - X^e"™) = 

One can see that the part of this equation corresponding to the symmetric nonlinear 
component of the Hamilton function depends only on {Xk,i} and the part corresponding 
to the asymmetric component can be expressed in {Xk,i} through (A3). Then we 
integrate this equation with respect to r 2 and obtain equation (11). 
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